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Abstract 

Thermodynamical properties of hot and dense nuclear matter are 
analyzed and compared for different equations of state (EoS). It is ar- 
gued that the softest point of the equation of state and the strangeness 
separation on the phase boundary can manifest themselves in observ- 
ables. The influence of the EoS and the order of the phase transition on 
the expansion dynamics of nuclear matter and strangeness excitation 
function is analyzed. It is shown that bulk properties of strangeness 
production in A-A collisions depend only weakly on the particular 
form of the EoS. The predictions of different models are related with 
experimental data on strangeness production. 
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1 Introduction 



The quest for the deconfinement transition, the phase transition from a 
confined hadronic phase to a deconfined quark-gluon phase (the so called 
quark-gluon plasma, QGP), remains a major challenge in strong interaction 
physics |rj . Over the past two decades a lot of effort has gone into the explo- 
ration of this transition and its possible manifestations in relativistic heavy 
ion collisions, in neutron stars as well as in the early universe. Relativistic 
heavy ion collisions offers a unique opportunity to reach states with tem- 
peratures and energy densities exceeding the critical values, T c ~ 170 MeV 
and e c ~ 0.6 GeV/fm 3 , specific for the deconfinement phase transition [5J. 
Thus, it is likely that color degrees of freedom play an important role al- 
ready at SPS and RHIC energies [3j. Various signals for the formation of 
a quark-gluon plasma in such collisions have been discussed and probed in 
experiments |TJ E] . 

Enhanced production of strangeness relative to proton-proton and proton- 
nucleus collisions was one of the conjecture signals of the quark-gluon plasma 
formation in heavy ion collisions [5J. The original idea behind the strangeness 
enhancement is that strange and antistrange quarks are easily created in a 
quark-gluon plasma, while in the hadronic phase strangeness production is 
suppressed. The dominant reaction in the plasma is gg — > ss. Furthermore, 
since the strange quark mass is not larger than T c , one expects the strange 
degrees of freedom to equilibrate in the quark-gluon plasma. Although, a 
heavy ion collision at high energies is a highly non-equilibrium process, the 
hadron yields (including strange particles) measured in the energy range 
from SIS to RHIC El El EE EI] are remarkably well described in the ther- 
mal model assuming chemical equilibrium at freeze out. This indicates that 
collective effects play an important role in the production of strangeness. 
On the other hand, elaborate microscopic transport models do not provide 
a quantitative explanation of the excitation functions for strange particles in 
this energy range. In the hadron string dynamics model JT] one finds a too 
small K + /tt + ratio around AGS energies, while in RQMD j^j the yield is 
overestimated at SIS and too small at SPS energies. 

The aim of this paper is to explore global effects of strangeness pro- 
duction in hot and dense nuclear matter within a collective approach. Our 
starting point is an equation of state (EoS) with a deconfinement phase tran- 
sition. Since strangeness is conserved at the time scales relevant for heavy 
ion collisions, a strangeness chemical potential is introduced. We examine 
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various phenomenological models for the equation of state, which differ in 
the order of the deconfmement phase transition: a first order transition (the 
two-phase bag model), a crossover-type transition (the statistical mixed- 
phase model) and no phase transition (pure hadronic models). The con- 
sequences of strangeness separation and softening of the equation of state 
are discussed. Furthermore, the manifestation of the order of the deconfme- 
ment phase transition in the expansion dynamics and the bulk strangeness 
production is studied. The predictions obtained with different equations of 
state are related with experimental excitation functions for relative strange 
particle abundances. 

2 Modelling the equation of state of strongly 
interacting matter 

The EoS of strongly interacting matter can in general be obtained by first 
principal calculations within lattice gauge theory ^3]. The thermodynam- 
ics and the order of the phase transition in QCD is rather well established 
for two and three light quark flavour in lattice calculations. However, the 
physically relevant situation of two light (u,d) and a heavy (s) quark is still 
not well described within lattice approach. In particular, the existence of a 
phase transition and its order in 2+1 flavour QCD is not yet known. In addi- 
tion most of the lattice calculations are performed for vanishing net baryon 
number density. Only recently, first results on the EoS with non zero baryon 
chemical potential have been obtained on the lattice [Tl]. However, these 
studies have so far been performed with large quark masses which distort 
the physical EoS. Thus, lattice results can still not be used directly in phys- 
ical applications. 

Lacking lattice QCD results for the EoS at finite baryon density % with 
physically relevant values of the quark masses, a common approach is to 
construct a phenomenological equation of state for strongly interacting mat- 
ter. This EoS should be constrained by existing lattice results and should 
also reproduce the two-phase structure of QCD. Here we construct different 
models for QCD thermodynamics and study their physical implications with 
particular emphasis on strangeness production and evolution in heavy ion 
collisions. 

A recent analysis of the lattice EoS shows that in the low tem- 
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perature phase, hadrons and resonances are the relevant degrees of freedom. 
The hadron resonance gas, with a modified mass spectrum to account for the 
unphysical values of the quark masses used in the lattice calculations, was 
shown to reproduce the bulk thermodynamic properties of QCD, obtained 
on the lattice with different numbers of quark flavors as well as at finite and 
vanishing net baryon density |To"l [TO] . 

Lattice calculations show that, at very large temperature the thermody- 
namical observables approach the Stefan-Boltzmann limit of an ideal gas of 
quarks and gluons, both at finite as well as vanishing net baryon density. 
The remaining ~ 20% discrepancy at T > 2Tq is understood by system- 
atic contributions in self-consistent implementation of quasiparticle masses 
in the HTL-reassumed perturbative QCD [T7j. To describe the thermody- 
namics near the phase transition additional model assumptions are necessary 

pang. 

From the above discussion it is clear that the most straightforward model 
for the EoS would be a non-interacting hadron resonance gas in the low 
temperature phase and ideal quark gluon-plasma in the color deconfined 
phase. These phases are matched at the phase transition boundary by means 
of the Gibbs phase equilibrium conditions. By construction, this approach 
yields a first order phase transition. Such an EoS with strange degrees of 
freedom has frequently been used in the literature [20J |2U 122 1231 12H I2H] and 
is also a standard input in hydrodynamic simulations of heavy ion collisions 
|2*m |27j. However, in order to obtain a reasonable phase diagram one has 
to include short-range repulsive interactions between hadronic constituents. 
In general this can be realized by introducing short-range repulsion in a 
thermodynamically consistent approach [2E1 EH USB EH] ■ 

We note that according to Gibbs phase rule [32], the number of thermo- 
dynamic degrees of freedom that may be varied without destroying the equi- 
librium of a mixture of r phases, with n c conserved charges is M = n c + 2 — r. 
For the hadron-quark deconfinement transition under consideration r = 2. 
If the baryon number is the only conserved quantity, n c = 1 and M = 1. 
Thus, the phase boundary is one-dimensional, i.e. a line. The Maxwell 
construction for a first order phase transition corresponds just to this case 
r = 2 and n c — 1. When both the baryon number and strangeness are 
conserved (n c = 2), one has M = 2 and therefore the phase boundary is in 
general a surface. In such a system, a standard Maxwell construction is not 
possible 

To account for the uncertainties in the order of the phase transition in 2+1 
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flavour QCD and also for the deviation of the equation of state from an ideal 
gas near the deconfinement transition we employ the EoS of the mixed phase 
model [2E1 12H] • In this model it is assumed that unbound quarks and gluons 
may coexist with hadrons forming a homogeneous mixture. This model is 
thermodynamically consistent and reproduces the lattice EoS obtained in 
the pure gauge theory as well as in two flavour QCD. Furthermore, the order 
of the phase transitions in the mixed phase model depends on the strength 
of the interaction between the phases. In this approach we can study the 
importance of the order of the phase transition on strangeness production 
and on the evolution of heavy ion collisions. 

In the following we discuss first the basic thermodynamical properties of 
these different models of the EoS and indicate relevant differences in their 
predictions. 



2.1 Two— phase bag model 



In the two-phase (2P) model 3 IJ. the deconfinement phase transition is de- 
termined by matching the EoS of a relativistic gas of hadrons and resonances, 
with repulsive interactions at short distances, to that of an ideal gas of quarks 
and gluons. The change in vacuum energy in the plasma phase is param- 
eterized by a bag constant B. We work in the grand canonical ensemble 
and account for all hadrons with mass rrij < 1.6 GeV, including the strange 
particles and resonances with strangeness Sj = ±1,±2,±3. The density of 
particle species j is then 



n 



(T, ft) = Uj(T, n B , /i s ) = v nf(T, fi B , fig) 



v 9 3 
2tt 2 



/ dk k 2 fj(k,T,fi B ,fi S ) , (1) 
Jo 



where 



fj{k,T, fiBifis) 



k 2 + m 2 - bjfiB - Sjfi S 
1 ■''!> I ~ j, I ± - ! - 



(2) 



is the momentum distribution function for fermions (+) and bosons (— ) while 
gj is the spin-isospin degeneracy factor. The chemical potential fij is related 
to the baryon (fj, B ) and strangeness (fig) chemical potentials 



fij = I), flu + Sj fi S 



(3) 
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where bj and Sj are the baryon number and strangeness of the particle. The 
quantity n* d corresponds to the number density of an ideal point-like hadron 
gas (IdHG). The factor 

v = v(T, n B , fx s ) = 1/[1 + X] v oj nf(T, fi B , Us)] (4) 

j 

reduces the volume available for hadrons due to their short range repulsion 
determined by the eigenvolume Voj = (l/2)(47r/3)(2r j) 3 We choose 

the effective interaction radius tqj ~ 0.5fm for all hadrons. Following (fl]). 
the baryon density and strangeness in the hadronic phase can be expressed 
as 

n B = J2 b i n j( T ^B,^s) , (5) 
jeh 

n s = J2 s i nj{T,fi B ,Hs) (6) 

j£h 

where the sum is taken over all hadrons and resonances. Similarly, the energy 
density of species j is given by 

£j(T, iab, Us) = v ef(T,HB,Hs) 

= m [°° dk k 2 Jk* + m) fj(k, T, hb, lis) ■ (7) 

Z7T Z JO 

In early studies j3H EH] , the excluded volume correction v was imple- 
mented in the same way for all thermodynamic quantities of the hadron gas, 
including the pressure 

P H (T,ii B ,Hs) =^2Pj(T,iiB,lis) , (8) 

jeh 

where the partial pressures are given by 

Pj(T,fi B ,fts) = v pf(T,ti B ,ns) 

= 6^1 dk ^/WT^ fj{k,T, ^ B ^ s) ' (9) 

However, this expansion for the pressure is not thermodynamically consistent 
with the charge (p3EJ) as well as the energy density 0. In Ref. [HE] & was 
shown that, it is possible to account for a thermodynamically consistent 
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implementation of the excluded volume corrections. In this approach the 
pressure is given by that of an ideal gas with modified chemical potentials 

p H (T^ B ^ s ) = J2pf(T,fI J ) (10) 

j eh 

where 

fij = bj [i B + Sj us - v 0j p H (T, n B , us) ■ (11) 

The remaining thermodynamic quantities are obtained with the excluded 
volume correction given above by taking the corresponding derivatives of the 
pressure. Thus, in this approach all fundamental thermodynamic relations 
are fulfilled We shall refer to Eqs. (|8I9|) and Eqs. (|10lll|) as two-phase 
thermodynamically inconsistent (2PIN) and consistent (2PC) model, re- 
spectively. Note that such an equation of state may violate causality at high 
densities, because an extended rigid body is incompatible with the basic 
principles of relativity. 

The QGP phase is described as a gas of non-interacting point-like quarks, 
antiquarks and gluons. The non-perturbative effects associated with confine- 
ment are described by the constant vacuum energy B. The pressure in the 
plasma is then given by 

p Q (T,fi B ^s) =p g (T) + J2P?(T,VB,fis) - B , (12) 
where the gluon 

P 9 (T) = (g g = 16) (13) 

and the quark pressure is obtained from Eq.fjHJ) for u,d,s quarks and anti- 
quarks. We use the quark masses m u = ma = 5MeV and m s = 150MeV 
and the bag constant B = (235MeV) 4 which yields a transition temperature 
T c 160MeV in agreement with lattice calculations at ub = JH]- The 
energy density of the plasma phase 

e Q (T, Li B , (is) = e g (T) + £ e) d {T, fi B , p s ) + B , (14) 
where the gluon contribution is given by 

e g (T) = 3 Pg (T) = ^T 4 (15) 
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and that of quark species j is obtained from Eq.(J7J) with v = 1. Analogously 
to Eqs.(j5J) and © the densities of the conserved charges in the QGP phase 
are : 

A = ^bjnfiT^B^s) , (16) 
n Q s = J2sjnf(T,tiB,Hs) ■ (17) 

The equilibrium between the plasma and the hadronic phase is determined 
by the Gibbs conditions for thermal (T® = T H ), mechanical {jft = p H ) and 
chemical {fi% = fig, /xf = fig) equilibrium. At a given temperature T and 
baryon chemical potential fi B the strange chemical potential fis is obtained 
by requiring that the net strangeness of the total system vanishes. Thus, for 
the total baryon density n B the phase equilibrium requires that: 

P H (T, Hb, fJ's) = P Q (T, [Lb, fJ's) , (18) 
n B = a n%{T,iJ lB ,ii s ) + (1 - a) n B (T,/j s ,/i S ) , (19) 
= a (T,fj. B ,fi s ) + (1 - a) nf (T,fj. B ,fi s ) , (20) 

where a = Vq/V is the fraction of the volume occupied by the plasma phase. 
The boundaries of the coexistence region are found by putting a = (the 
hadron phase boundary) and a = 1 (the plasma boundary). 

As mentioned above, the Maxwell construction is not appropriate in a 
system where both baryon number and strangeness are conserved. To il- 
lustrate this, we first analyze an approximate form of the equation (|17|) for 
strangeness conservation. We retain only the main terms and drop those 
with \sj\ > 1: 

a (n s - n s ) = (1 - a) (n K + n A + n t - n R - n A - n s ) . (21) 
In the Boltzmann approximation the densities may be computed analytically 

nf - nf = 9j (^) (^f) 2 K 2 (^-) exp(|) = 9j (ll) Wj exp(|) , (22) 
and the strangeness chemical potential is obtained as 

Vs = -\n-— j L- — — — U tf. (23) 



— in 

2 3aW s + v(l - a) [W K e t# + (W A + 3l¥ s )e 
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It is seen that at the plasma boundary (a = 1) fig = /Ab/3 while fig ^ 
at the hadron boundary (a = 0). This implies that not only fig but 
also fiB and the pressure changes along isotherms in the coexistence region. 
Hence, the standard Maxwell construction, which interpolates the densities 
linearly between the pure phases, is not adequate. The equations for phase 
equilibrium (|18H2(J|) must be solved to obtain fig and at every point in 
the coexistence region. 

When two phases coexist, the system is in general not homogeneous as 
the phases occupy a separate domains in space. We do not explicitly account 
for such domains structure nor for a possible surface energy contribution 
to the equation of state. The only consequence of the phase separation in 
this calculations is that the interactions between particles in the plasma and 
hadronic phase are excluded. This is different in the statistical mixed phase 
model discussed in the next section. 

The solution of the Gibbs conditions (|18H20p is shown in Fig. ^ for the 
plasma and hadron phase pressure versus fi% at fixed T = 80MeV and 
fig = yU_e/3. The crossing of the quark and hadronic pressure corresponds to 
the transition point at the plasma boundary. In this special case the condi- 
tion fig = guarantees strangeness neutrality. In general, however, for 
a 7^ 1, fig must be chosen such that the strangeness of the total system of 
quarks and hadrons vanishes. This requires an iterative solution of the equa- 
tions (jl8H2(Jj) . Away from the transition point, the system is in the phase 
with higher pressure p (lower free energy). Fig. Q also shows that there is 
no deconfinement transition if the hadronic phase is described as a gas of 
point-like particles jSl]. The situation is not improved by including more 
resonances. On the contrary, the larger the set of hadronic resonances is, the 
higher is the pressure at a given baryon chemical potential. However, the 
inclusion of repulsive interactions between hadrons leads to a reduction of 
the hadron pressure p H at fixed baryon chemical potential. Consequently, a 
short-range repulsion between hadrons stabilizes the quark-gluon plasma at 
high densities. 

The resulting phase boundaries in the T-fi plane are shown in Fig. |21 The 
difference in fi B at the phase boundaries described by Eqs. (fTHH2T!|) is small 
while for the strange chemical potential fig it is more noticeable. It is natural 
to expect that in the high temperature plasma fig ~ /is/3. On the other 
hand, in the hadronic phase and at low temperatures, where strangeness is 
carried mostly by kaons and A-hyperons, the strange chemical potential is 
roughly approximated by fig ~ 0.5 (fiB + m K — m \) ~ 550MeV. Both these 
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Figure 1: Pressure versus baryon chemical potential for fixed T = 80 MeV 
and for fi$ = The thin-line is the hadronic and the thick-line the 

quark phase in 2PC model. The dashed-dotted (1) line and dashed-line are 
an ideal gas model results without and with repulsion in the 2PIN model, 
respectively. The line (2) is obtained as line (1) but with fewer hadronic 
resonances. The line (3) is calculated within a mean-field approximation of 
the Zimanyi model |37j (see text). 



expectations are in agreement with our numerical results. Nevertheless, also 
in the high temperature hadronic phase the strange chemical potential ex- 
hibits an approximately linear dependence on the baryon chemical potential. 

In Figs. El and El the resulting phase diagrams are shown in the T—fig 
and T-/ig as well as T-i%b planes. The role of thermodynamical consistency 
is particularly evident in the T-rig plane. As seen in Fig. EJ the baryon 
density «b at the plasma boundary is increased while it is slightly decreased 
at the hadron side in the 2PC model as compared with the 2PIN approach. 
Consequently, the range of the coexistence region grows from ~ (4 -f- 7.5)n 
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Figure 2: The phase diagram in the T-/zb (marked by B) and in the T-fi$ 
(marked by S) plane for the 2PC and 2PIN models. The plasma and hadron 
boundaries are shown by full and dashed-lines, respectively. The dotted- 
lines are the approximate results obtained with fi S = and with /i S from 
Eq.(J2i. 
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to ~ (3.5 + 10)n . 



2P-consistent model 



T=U0 MeV 




n B /n n B /n 

Figure 4: Dependence of different thermodynamical quantities on the baryon 
density within 2PC model. The results are shown for two different tempera- 
tures. The hadron and plasma phase boundaries are shown by dotted-lines. 

Thermodynamical properties and the differences between the two-phase 
bag models are shown in Figs. 0] and El In both cases the baryon and strange 
chemical potentials are continuous when crossing the phase boundaries. This 
guarantees that the system is chemically stable. Demanding the conservation 
of strangeness in each phase separately [20] would results in a discontinuity 
in us ■ 111 contrast to the case with only one conserved charge, the chemical 
potentials are not necessary constant within the Gibbs coexistence region. 
Depending on the values of [i at the hadronic and plasma boundaries (see 
Fig. |2J), the chemical potentials (in particular \xs) can be either increasing or 
decreasing functions of ns- Although this change is not large, it influences 
the strangeness separation in the phase coexistence region. 

The energy density is seen in Figs. H] and El to be a monotonously in- 
creasing function of in both models. The pressure is also continuous 
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within 2PC model and is higher than in the 2PIN approach. In addition, 
in the latter model the pressure also suffers a jump at the boundary of the 
hadronic phase, which increases with decreasing temperature. Such an EoS 
would lead to a mechanical instability of the hydrodynamic flow. As seen 
in Figs. 0] and the changes in pressure across the coexistence region are 
quite small. Consequently, the system expands very slowly. This is a specific 
feature expected for the systems with a first order phase transition. 



T ' 1 " 1 " 1 " — I — I ' 1 " 1 

2P-inconsistent model 




2 l 6 2 

n B /n n B /n 



Figure 5: The same as in Fig. H]but for the 2PIN model. 

We stress that there are at least two problems, which show up when 
the EoS discussed above is employed in hydrodynamic calculations. First, 
as shown in EE] causality is violated at densities ub > 3.5n . Second, 
the ideal gas model with an excluded volume correction does not repro- 
duce the saturation properties of nuclear matter. An attempt to combine 
the excluded volume correction with a mean field treatment of the hadronic 
interactions resulted in an incompressibility parameter which is too large 
K > 550 MeV jH]. 
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2.2 Statistical mixed phase model 

The mixed phase (MP) model [2H1 EH EH] is a phenomenological model 
of the EoS with a deconfinement phase transition of QCD which shows a 
satisfactory agreement with the lattice data. The underlying assumption of 
the MP model is that unbound quarks and gluons may coexist with hadrons 
forming a spatially homogeneous quark/gluon-hadron phase which we call a 
generalized Gibbs mixed phase. Since the mean distance between hadrons 
and quarks/gluons in the mixed phase may be of the same order as that 
between hadrons, the interactions between all these constituents (unbound 
quarks/gluons and hadrons) play an important role. The strength of this 
interactions defines the order of the phase transition. 

To find the free energy within the MP model |25| , the following ef- 
fective Hamiltonian, expressed in terms of quasiparticles interacting with a 
density-dependent mean field, is used : 

# = ££/dr^+(r,a) (/ : VH^+^(p)) Uv,a)-C(p)V .(24) 

Here ipi(r, a) denotes a field operator for the quasiparticle species i charac- 
terized by the mass (the current masses for quarks and gluons and the 
free hadron masses are used here). The index a accounts for spin, isospin 
and color degrees of freedom. Furthermore, U{ is the mean field acting on 
particles of type i, C(p) is a potential energy term, which is needed to avoid 
double counting of the interaction, and V is the volume of the system. 

By requiring thermodynamical consistency j2Hl 123 EDI lEJj one finds con- 
straints on the parameters in the Hamiltonian. The constraints follow from |2EJ 

EH] 

<§>=0, <^> = 0, (25, 

where (. . .) denotes the statistical average. For the Hamiltonian (|24j). these 
conditions reduce to 

dUi dC dU t DC 

^ p *dp- - dp- = ' £>ar - &T = °' (26) 

which, as shown in [2H1 12H] , imply that Ui(p) and C(p) do not explicitly 
depend on temperature. 
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We model color confinement by assuming the following density depen- 
dence for the mean-field potential of quarks and gluons 

U q (p) = U g (p) = ^ ; 7 >0 (27) 

where 

P = P q + Pg + Y.Pj = Pi + P9 + Y. u 3 n i ( 28 ) 

j 3 

is the total number density of quarks and gluons in the local rest frame and p q 
and p g are the number densities of unbound (deconfined) quarks and gluons 
(p p i = p q + p g ), while rij is the number density of hadrons of type j having Vj 
number of valence quarks inside. The presence of the total number density p 
in ()27|) implies interactions between all components of the generalized Gibbs 
mixed phase. The potential (J27|) exhibits two important limits of QCD. 
For p — > 0, the interaction potential approaches infinity, i.e. an infinite 
amount of energy is necessary to create an isolated quark or gluon. This 
obviously simulates confinement of colored objects. In the opposite limit 
of large energy density, p — > oo, we have U g — > which is consistent with 
asymptotic freedom. 

In the description of the hadron components, the MP model accounts not 
only for hadron-hadron but also for quark/gluon-hadron interactions. The 
mean field acting on the hadron species j in the MP model has two terms : 

U j = uf + uf l) . (29) 

In the limit where there are no unbounded quarks and gluons, Uj = 0, i.e., 

Uj = Uj h \ This happens at low densities, where colored degrees of freedom 
are confined in hadrons. 

Due to the constraints the second term in Eq. (j29|) my be written as 
EEj: 

U?> = ^(l-(l-M- 7 ) , (30) 

where w p i = p p i/p is the fraction of quark-gluon plasma in the mixed phase l . 
Thus, if U q and U g are known, the thermodynamic consistency conditions f!26D 
allow us to unambiguously determine the correction term C(p) in Eq. ([23$ . 

1 We note that the resulting hadron single particle potential becomes very attractive 
near the plasma phase boundary, and even diverges in the plasma phase. We do not believe 
that this behaviour is physical. Nevertheless, we adopt this convenient prescription, since 
it does not have a noticeable effect on the thermodynamics nor on the hydrodynamic flow. 
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The hadronic potential Uj is described by a non-linear mean-field mo- 
del §7\ 

u j = 9r,j + g a ,j ^2(2/) , (31) 

where g r j > and g a j < are repulsive and attractive coupling constants, 
respectively. 

Thermodynamic consistency implies that the functions tpi(x) and ^2(2/ ) 
depend only on particle densities. In Ref. |37] these functions are chosen 
such that 

b 1 (p 1 = x, — 6i(y 2 + hvl) = y (32) 

where 

X = J2 9r,i Pu,, V = J2 9a,i p Vi ■ 

and b\ and b 2 are free parameters. In j37j considering a mixture of nucle- 
ons and A's the model parameters were fixed such that to reproduce the 
saturation properties of nuclear matter and the ratio of the A to nucleon 
coupling constants. We generalize this approach by including all hadrons in 
our model and assuming that the coupling constants scale with the number 
of constituent quarks : 

U j k) = v i (^i0° ~ Pvi) + <?2{P~ Ppl) ) , (33) 
where (p\ and (p 2 satisfy the equations 

ci<fi=p-p p i, -c 2 <f 2 -c 3 ifl = p-p p i (34) 
with p — ppi = Y, u - u jPj- The parameters in Eq. are given by |2H] 

bi _ 61 b x b 2 

and are fixed by requiring that the properties of the ground state (T = 0, 
tib = n o ~ 0.17 fm~ 3 ) of nuclear matter are reproduced: i.e. a binding 
energy per nucleon of —16 MeV, incompressibility of 210 MeV and vanishing 
pressure. 

We also addressed to the extension of the Zimanyi model [37J as interact- 
ing hadron gas (InHG) model with no phase transition. The /^-dependence 
of the pressure in this model is illustrated in Fig. ^ 
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The thermodynamics in the MP model is obtained from the partition 
function in the standard way. The baryon and strange chemical potentials 
are fixed by the baryon number and strangeness conservation, 

n B (T,fi) = Y, h J "./^ '/'•/'/) • (35) 
n s {T,fi)=0 = s i n j( T > Vj) , (36) 

where the sum is taken over all quarks, gluons and hadrons. The same set 
of hadrons and resonances is used here as in the previous models. 
As an example we quote an expression for particle number density, 



nj(T,fij) = ^ J dkk 2 

In the energy density and pressure there are further terms, originating from 
the interactions. These terms, are obtained |28[l2"D] by solving the consistency 
conditions (J2HJ). 

2.3 Thermodynamics in the mixed phase model 

The parameters of the MP model are fixed [2E1 HHj by requiring that the lat- 
tice results for thermodynamical quantities in the pure gauge and two flavor 
QCD are reproduced. In this study the quark-gluon interaction parameters 
7 = 0.62 and ^.Vw-t-i) = 250 MeV were obtained. The same parameters are 
assumed to be valid also for the interactions with strange quarks. 

In Fig. El we compare the thermodynamical properties of the MP and 
2PC models. The energy density and pressure in the 2PC model shows the 
typical behavior for a system with a first order phase transition: an abrupt 
change in the energy density at T = T c and smooth change in the pressure. 
In the MP model on the other hand, both e and p vary continuously with 
temperature. The transition temperature in the MP model, T c « 160 MeV, 
is defined by the maximum of the heat capacity (see insert in Fig. [HJ). 

In Fig. [7|the ratio p/e is shown for different values of the total baryon 
number density in three models under consideration. A common feature 
of these models is that for finite baryon densities they all exhibit a clear 
threshold behavior. The threshold is shifted to higher e with increasing 
However, in contrast to the ideal hadron gas, both the MP model and 2P 



exp 



'k 2 + m 2 — bjfi B 



± 1 



.(37) 
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Figure 6: Temperature dependence of the energy density and pressure at 
vanishing total baryon density. Full and dashed-lines are the MP and 2PC 
model results, respectively. The insert figure shows the reduced heat capacity. 



model have the softest point in the EoS, i.e., a minimum [26 of the function 
p(e)/e. A particular feature of the MP model is that even for % = the 
softest point is not very pronounced and located at a relatively low energy 
density: esp ~ 0.45 GeV/fm 3 . This is consistent with lattice result jlUj . 
In the MP model, the softest point is gradually washed out with increasing 
baryon density and vanishes completely for ns > 0.5 uq. This is, however, 
not the case in the 2P models, where one finds a pronounced softest point at 
large energy density e$p ~ 1.5 GeV/fm 3 , which depends only weakly on the 
baryon density % see Fig. |7J). Finally, in the InHG model as well as in the 
relativistic ideal hadron gas there is obviously no softest point in the EoS. 

The differences in the thermodynamical properties of the above models 
will be also reflected in the expansion dynamics of a thermal fireball created in 
heavy ion collisions. The effect of these differences on strangeness production 
and evolution will be explored in the following sections. 
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Figure 7: The ratio of the pressure (p) to the energy density (e) as a function 
of e. The results are for different values of a total baryon density (ub) and 
for three models of the EoS. 



3 Strangeness production 

3.1 Strangeness content in equilibrium 

The conservation of strangeness in the coexistence region of quarks and 
hadrons implies that the total number of strange and antistrange quarks 
are equal. However, the s—s content in the individual phases may differ from 
zero. The strangeness content of the quarks in the mixed or plasma phase 
is characterized by two ratios : p s / Ps and D s = (p s + ps)/p P i (see Eq. (j2Hj) ). 
The second ratio gives the strangeness fraction in the plasma. 

In Fig. IHlthe ratio p s j p s is shown as a function of ps for a fixed plasma 
fraction a. For a ~ 1 the ratio p s j p s ~ 1 for almost all values of ps- However, 
if a << 1, that is when the volume in the mixed phase is mostly occupied by 
hadrons, the separation of strange and antistrange quarks is clearly seen in 
Fig. |H| This is mainly because the hadronic component of the mixed phase is 
dominated by the kaons, while the hyperons are suppressed due to their large 
masses. This strangeness excess through kaons is compensated by creation of 
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s-quarks in the plasma. The results in Fig. |S] are in qualitative agreement 
with Ref. j22] where the 2PIN model without higher mass resonances was 
employed. The contribution of higher mass resonances results in an increase 
of p s /ps for (/x b /3)h « 400 - 500 MeV. 
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Figure 8: Ratio of strange to antistrange quark densities in a quark-gluon 
plasma calculated along the hadronic boundary. The results are for PC 
model calculated with different values of the volume fraction (a) occupied 
by a quark-gluon plasma. 

In Fig. El the strangeness composition in an equilibrium system is com- 
pared for two different models. In the bag model EoS and at high temperature 
(T ~ 140 MeV) the p s j p s ratio decreases when the baryon density inside the 
Gibbs mixed phase approaches the plasma boundary. However, for the mod- 
erate temperatures (T ~ 80 MeV), the ratio p s j ' p s < 1 and it increases with 
ub- The above behavior is a direct implication of the simultaneous conser- 
vation of strangeness and the baryon number. If these conservation laws are 
decoupled [22], then this behavior at low temperatures is not seen. 

In the MP model the p s j p s > 1 for all values of the baryon density. For a 
fixed temperature the p s j ' p s ratio is seen in Fig. EJto gradually decrease with 
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Figure 9: The p s j p s ratio for quark component and strangeness fraction (D s ) 
for unbound quarks as a function of baryon density. The results are shown 
for different temperatures and for two EoS. The plasma boundary is marked 
by arrows. Note the factor 1/5 in the MP model at T = 80 MeV. 



increasing density. Its values are noticeably higher than in the 2P model. In 
both models, however, the strangeness separation effect is stronger when the 
system is closer to the hadronic boundary, i.e. where there is small admixture 
of quarks. For the 2P model this corresponds to the existence of a small blob 
of plasma while in the MP model a homogeneous admixture of unbound 
quarks and gluons with small concentration. 

Above the hadronic phase boundary, the n^-dependence of D s in the 
2PC model is similar to that in the MP model. The strangeness fraction in 
the MP model is the largest below the hadronic boundary and maximal in 
baryon free matter. In Fig. El we note a jump in Ds which corresponds to 
a jump in strange particle multiplicity when crossing the phase boundary; a 
similar jump is observed in the baryon number. 

From the above discussion, it is clear, that the strangeness content and 
its distribution in the transition region from the quark-gluon plasma to the 
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hadronic phase is strongly model dependent. It is effected by the order of the 
phase transition and the strength and the form of the interactions between 
constituents. These differences are particularly evident at moderate values 
of the temperature and baryon density. This is just the region which is 
traversed by an expanding system created in heavy ion collisions on its way 
towards the chemical freeze-out. Thus, one could expect that the order of the 
phase transition and particular strangeness dynamics could manifest itself in 
observables in heavy ion collisions. 

3.2 Strangeness evolution in expansion dynamics 

To study the possible influence of the EoS on observables in heavy ion colli- 
sions, we have to describe the space-time evolution of a thermal medium that 
is created in the initial state. This is conveniently done within a hydrody- 
namical model. The EoS is an input for constructing the energy — momentum 
tensor, which is needed in the hydrodynamical equations. 

To solve the hydrodynamical equations for a given experimental set-up 
one needs to specify the initial conditions. The initial volume, the entropy 
and the baryon number densities in the collisions are modelled within QGSM 
transport code jH]. The predictions of this model are consistent with the 
results obtained within the RQMD and UrQMD transport codes. 

We assume that, in the center of mass frame, the initial state is a cylinder 
of radius R = 5 fm and Lorentz contracted length L = 2R/ r y c _ m _. This initial 
state corresponds to the time when the centers of the colliding nuclei just 
have passed the point of full overlap. 2 We neglect the transverse expansion 
and assume that the hydrodynamical evolution of the fireball is described by 
a one-dimensional isentropic expansion of the scaling type in the longitudi- 
nal direction. In this approximation the entropy and baryon density decrease 
inversely proportional to the expansion time. The values of all other ther- 
modynamic quantities are obtained from the EoS at each temporal step (see, 
for example [42J). 

In Fig. El we show the fireball evolution trajectories for central Au-Au 
collisions in the T— plane for different collision energies and for different 
EoS. The chemical freeze-out parameters obtained [HI El EH El] within the 
statistical model at different collision energies are also shown in this figure. 

2 A detailed description of the procedure to fix the initial conditions in heavy ion colli- 
sions can be found in fM 1351 1ST] . 
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Clearly, the chemical freeze-out parameters from SIS up to RHIC are well de- 
scribed by the universal condition of fixed energy/particle, (E had ) / (N had ) ~ 
lGeV@3IE5|. 
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Figure 10: A compilation of the chemical freeze-out parameters from Refs. 
[3 EH Hlj obtained with the hadron resonance gas partition function at dif- 
ferent beam-energies (filled dots, squares and triangles). The smooth dashed 
curve is the universal freeze-out curve of fixed (E had ) / (N had ) ~ 1 GeV from 
Ref. |43j . Also shown are dynamical trajectories for central Au-Au collisions 
calculated within different models (the interacting hadron gas model (InHG), 
the mixed-phase model (MP), the ideal hadron gas model (IdHG) and the 
thermodynamically consistent two-phase model (2PC)). The empty circles 
near the end of each trajectory correspond to freeze-out condition of fixed 
energy density, £/ — 0.135 GeV/ fm 3 . 

The dynamical trajectories show a strong dependence on the properties 
of the EoS. In the MP model there is a turning point seen in all trajectories, 
i.e. the point where dT/dfis changes sign. The existence of such a point is 
a general feature of the MP model and is directly related to the appearance 
of two limiting regimes: 

(i) At high temperatures and in the ultra-relativistic limit, m q — > 0, 
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the thermodynamic potential = — Vp can be obtained analytically from 
Eqs.© an d (HI 



VL=-V ( ai T 4 + a 2 T 2 fi% + a 3 /4) . (38) 
The entropy per baryon 

s on/or 2a ± + fl 2 yjr) 

- B -dmr B ~ a2 { ^ )+as ( /^ )3 ' (39) 

is conserved along trajectories defined by fi B /T = const. Thus, in the high 
temperature limit, an isentropic expansion is characterized by a linear rela- 
tion between T and fi B . 

(ii) At intermediate temperatures, the system can be approximated by a 
Boltzmann gas (|22j) of a non-relativistic nucleons. In this case the entropy 
in the dilute gas approximation is given by 

S = -J^- 3 Jd 3 P [/ln/ + (l-/)ln(l-/)] 



N 



B 



Jd 3 p f\nf 
JdPpf 



(40) 



with the distribution function / = exp[(/i B — m B — p 2 /2ms)/T]. In this 
temperature range, conservation of the entropy per baryon implies that 

^ = h - + ^^ = const. (41) 
n B 2 T K J 

Thus, for intermediate temperatures, we again find a linear relation between 
T and \i B but with a negative slope. The different behavior of \i B (T) at high 
and intermediate temperatures, implies that there is a turning point in the 
fireball expansion trajectories, as seen in Fig. 

The dynamical trajectories calculated in the MP model pass quite close to 
the phenomenological freeze-out points. For all collision energies, the turning 
point is located on the universal freeze-out curve of fixed energy/particle. 
This fact has been noticed already in Ref. jlH] for the MP model with 
two light quarks. The contribution of strange quarks and the requirement 
of strangeness conservation modifies the dynamical expansion path of the 
fireball. This is particularly evident for Ei a b < 10 AGeV where neglecting 
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the strange quarks gives rise to a visible shift of the turning point towards 
smaller p B . 

In the parameter range below the phenomenological freeze-out curve the 
expansion path in the MP, IdHG as well as in MP model are quite similar. 
In the InHG model, however, there is a small shift toward larger values of 
fis- This agreement indicates that in the final stage the expansion path 
depends only weakly on details of the equation of state. The dynamical path 
is, to a large extent, determined by the entropy/baryon and strangeness 
conservation which in the hadronic phase puts strong constraints on the 
particle composition of the fireball. In this case the space time evolution and 
thermodynamics is governed by a gas of weakly interacting resonances, the 
effective degrees of freedom in the low temperature phase of QCD. This may 
be the reason behind the success of the non-interacting hadron resonance 
gas in the description of bulk observables in heavy ion collisions. 

The differences between the various equations of state in the evolution of 
the thermal fireball is clearly visible above the freeze-out curve. In contrast 
to the MP model, the IdHG turning points do not correlate with the freeze- 
out curve. There is also no softest point in the InHG and IdHG model. The 
dynamical trajectories within the 2P bag-models exhibit a characteristic re- 
heating regime in the phase coexistence region. For this model, the expansion 
trajectory closely follows the phase boundary in this regime, as shown in 
Fig. |21 At SPS energies and above, the hadronic end of the intermediate 
coexistence region in the T—pis plane (the so-called "hottest hadronic point") 
is close to the phenomenological chemical freeze-out point. At lower energies 
there is no such correlation for the 2P models. For Ei a b < 10 AGeV the initial 
state is in the phase coexistence region. 

The question of strangeness separation in heavy ion collisions addressed in 
Ref. |221 f° r static system can be reanalyzed in our approach for dynamically 
evolving fireball. The results are shown in Fig. ^2 for Au-Au collisions at 
different bombarding energies within 2P and MP models. In both models 
Ps/ps > 1, since there is no chance for the system to pass through a high 
density baryonic state where p s could be less then unity. 

In the 2P model we find that strangeness is separated to a less degree 
at the exit point from the phase coexistence region than found in j22]- On 
the other hand in the MP model the system evolves much longer and conse- 
quently a higher degree of strangeness separation is obtained. This effect is 
stronger at E lab = 10 AGeV than at 160 AGeV. 

So far the differences between various models for the expansion dynamics 
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Figure 11: Time evolution of the ratio of strange to anti-strange quark 
densities in the hadronic component. The MP and 2P model results are for 
central Au-Au collisions at different beam-energies. 



were discussed on the level of global thermodynamical quantities. It is of 
particular interest to explore physical observables that are directly measured 
in heavy ion collisions. In the following we consider strange particle mul- 
tiplicity ratios to discuss the influence of the equation of state on particle 
yields. The predictions of the different models will be compared at thermal 
freeze-out where the particle momentum distributions are frozen. We assume 
a shock-like freeze-out j37] where energy, the total baryonic and strangeness 
charges are conserved. The thermal freeze-out conditions are assumed to be 
determined by the fixed energy density ej ~ 0.9n mjv = 0.135 GeV/fm 3 . 
Below this energy density the system consists of a free streaming gas of non 
interacting particles. The thermal freeze-out points are shown in Fig. [TUlbv 
empty circles on each trajectory for all models and for all collision energies. 

The excitation function of the relative yields of K + mesons calculated 
within the MP model is shown in Fig. as triangles. For reference, we also 
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show in this figure the 47r-integrated data for K + /tc + ratio obtained in heavy 
ion collisions at different beam energies. The shape of the kaon excitation 
function in the MP model is similar to that seen in the data. However, the 
absolute values are overestimated, especially for the low collision energies. 
We have to stress, however, that the models discussed here are still not quite 
suitable to be compared with data. First, the conservation of electric charge 
was not taken into account. The isospin asymmetry is particularly relevant 
at low collision energies (below AGS) where it can change the charge par- 
ticle multiplicity ratios by up to 20%. Second, the hydrodynamical model 
applied here describes a longitudinally expanding fireball. This is, to a large 
extend sufficient at RHIC energies, however, it may be not valid at AGS or 
SIS where transverse expansion cannot be neglected. Furthermore, only part 
of particle mass spectrum was included with the masses up to 1.6 GeV. At 
AGS and higher energies the contributions from heavier resonances increases 
the yields of lighter particles. Finally, the system may be out of chemical 
equilibrium at some stages during the evolution from chemical towards ther- 
mal equilibrium [3B]- Nevertheless, all these effects cannot account for the 
observed discrepancy by a factor of five between the MP model and data at 
low collision energies (Fig. [T2*|) . However, the differences may be due to the 
grand canonical (GC) treatment of the strangeness conservation used in the 
calculations. 

In the GC ensemble strangeness is conserved on the average and is con- 
trolled by the strange chemical potential. Within the statistical approach, 
the use of the grand canonical ensemble for particle production can be justi- 
fied only if the number of produced particles that carry a conserved charge 
is sufficiently large. In this case also event-averaged multiplicities can be 
treated in a grand canonical formulation. In this approach, the net value of 
a given charge (e.g. electric charge, baryon number, strangeness, charm, etc.) 
fluctuates from event to event. These fluctuations can be neglected (relative 
to the mean particle multiplicity) only if the particles carrying the charges in 
question are abundant. Here, the charge is indeed conserved on the average 
and a grand canonical treatment is adequate. However, in the opposite limit 
of low production yields (as is the case for strangeness production in low en- 
ergy heavy-ion collisions) the particle number fluctuation can be of the same 
order as the event averaged value. In this case charge conservation has to 
be implemented exactly in each event [HI ES] • In the statistical physics the 
exact conservation of quantum numbers requires a canonical (C) formulation 
of the partition function. 
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Figure 12: The K + /tt + ratio as a function of a beam-energy. The data 
points are from The lines are the MP model results obtained in the 

grand canonical as well as canonical formulation of strangeness conservation 
for different parameterizations of the volume parameter V c (see text). 



The grand canonical Z GC and canonical Zg partition functions are con- 
nected by a cluster decomposition in the fugacity parameter (A = exp(/i s /T)) 

s=oo 

Z GC (T, V, n B , A) = £ X s Z G (T, V, hb) . (42) 

s=— oo 

The relation (J42|) can be inverted and the canonical partition function with 
total strangeness S = is obtained from 

Z° =0 (T, V, pb) = -L \^ d<\> Z GC (T, V, fji B , A - e«) . (43) 

Z7T JO 

Neglecting the contributions from multistrange hyperons and assuming Boltz- 
mann statistics the density of kaons in the C ensemble is given by [201 03] 

n c - n B S ± (U) 
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where the argument of the Bessel function I s (x) is 



x = 2^SiS_i. (45) 

with 



B 



Here the particle density for hadron species j is given by Eq. (|22|) with 
— A*b^- The sum is taken over all particles and resonances carrying 
strangeness s. The volume V c is a model parameter which is interpreted 
as the strangeness correlation volume. 3 In the equilibrium statistical model 
a correlation volume V = V\ ~ 1.9nA part /2 was found to reproduce the 
experimental multiplicity ratios for all measured particle. In our dynamical 
approach, V c is assumed to be the initial volume of the collision fireball 
Vc = Vo(Ei a b), and thus is energy dependent [H2~] . 

From Eqs. (j44j) and (|2*2*|) it is clear that grand canonical and canonical 
results for the kaon density are related by the substitution [7j: 

exp( ^ /T) vlfc W) ■ (46) 

Thus, the main difference between C and GC results is contained in a reduc- 
tion of the fugacity parameter by the factor F = Ii(x) / 'Iq(x). In the limit of 
a large volume V c , i.e., x — > oo, the ratio F — > 1 and the GC and C results 
coincide. However, in the opposite limit, x — > 0, the factor F < 1 which 
leads to a suppression of the strange particle densities. The canonical sup- 
pression depends strongly on the temperature and the correlation volume. 
Both these parameters, in particular the temperature, are dependent on the 
collision energy. 

The sensitivity of the suppression factor on Ei a b and the correlation vol- 
ume is shown in Fig. The canonical suppression factor calculated with 
V c = V\ is increasing with collision energy and reaches its asymptotic value 
at Ei ab > 10 GeV. Obviously the magnitude of the suppression is strongly 
dependent on the correlation volume. This is shown in Fig. f° r three 
different parameterizations of V c . Particularly interesting is the behavior of 
F for V c = Vo(Ei a b). Due to the Lorentz contraction the initial volume is 



3 For a more detailed discussion of the interpretation and the role of this parameter 
see e.g. Ref. 
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Figure 13: The beam-energy dependence of the strangeness suppression fac- 
tor for central Au-Au collisions. This factor is calculated within MP model 
at freeze-out for different parameterizations of the correlation volume V c (see 
text). 



decreasing with increasing collision energy. At lower collision energies this 
decrease of the volume is compensated by an increase of temperature such 
that the suppression factor increases with E[ ab . However for E iab > 10 GeV 
there is only a moderate increase of freeze-out temperature that is not suf- 
ficient to overcome a decrease of V®. Consequently for Ei ab > 10 GeV the 
suppression factor starts to decrease with energy. 

The amount of canonical suppression at fixed Ei ab also depends strongly 
on the temperature which in turn is determined by the energy density at 
freeze-out. In the equilibrium analysis of particle production at SIS [7j the 
energy density at chemical and thermal freeze-out was a factor of three lower 
than the value used in the present dynamical study, e ~ 0.135 GeV/fm 3 . 
Consequently, for V c = Vq and for 1 < Ei ab < 2 GeV the canonical suppression 
found in Ref. [Zj was much stronger than that shown in Fig. ED We have 
not tuned the parameters to reproduce previous results. In low energy heavy 
ion collisions the expansion trajectories and the freeze-out parameters will 
change once the transverse expansion is taken into account. 

In fig. El we show the effect of the canonical suppression on the K + 
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excitation function calculated in the MP model with two different parame- 
terizations of the correlation volume: V c = Vq and V c = min(Vo,X^) where 
V2 = Vi/5. As expected, there is a noticeable decrease of K + yield due to the 
exact treatment of strangeness conservation. The suppression of strangeness 
at energies beyond AGS is entirely due to the energy-dependent Lorentz 
contraction of the initial correlation volume. In Fig. [Tlflthe results of a cal- 
culation also are presented where the choice of V c — min(Vo, V 2 ) is optimized 
to reproduce the K + /n + data. 




Figure 14: The ratios of 47r-integrated strange particle yields per pion yields 
for central Au-Au collision as a function of beam-energy. The compilation of 
experimental data is taken from OEH]- The calculated excitation functions 
are for different EoS with the canonical suppression factor. 



The above analysis of K + excitation function clearly shows that due to 
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associated strangeness production and the small production cross sections at 
low collision energies one has to implement exact strangeness conservation. 
In the following we will implement this concept in all models and discuss 
the predictions for strangeness production and energy dependence. In the 
calculations we use the correlation volume V c = min(Vo, V2). 

In Fig. El w e calculate relative excitation functions for different strange 
mesons and baryons for four hydrodynamical models. The most striking 
result seen in this figure is that all models yield very similar results for the 
strangeness excitation functions. This is particularly true for the production 
of K + /tt + and A/tt + where the results of all models besides InHG, are hardly 
distinguishable. Some differences are seen on the level of K~ excitation 
function which are mainly due to larger sensitivity of K~/tt~ ratio to the 
value of the temperature. It is interesting to note that all models show a 
maximum in the A/n excitation function for 10 < Ei a b < 30 GeV. Such a 
maximum is found also in equilibrium models [53] . 
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Figure 15: The Wroblewski ratio X s as a function of beam-energy for cen- 
tral Au-Au collisions. The contributions of mesons and baryons are shown 
separately. The points at AGS energies are from [44. 

The relative strangeness content of the produced particles in heavy ion 
collisions is characterized by the Wroblewski factor [HU EH] , 
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< uu > + < dd > 

where the quantities in angular brackets refer to the number of newly created 
quark-ant iquark pairs. The Wroblewski factor is shown in Fig. [T51for differ- 
ent collision energies. The separate contributions to \$ from strange mesons 
and baryons as well as its overall value is calculated within the MP, 2P and 
IdHG models. The results are compared with As obtained in an equilibrium 
model analysis of experimental data at AGS energies. There is a surprising 
agreement of all dynamical models on the relative strangeness content of the 
fireball at freeze-out. The results are also consistent with the equilibrium 
model [S^. However, the maximum spread of the Wroblewski factor seen in 
Fig. is broader than previously seen in the equilibrium canonical model 
|54| . In the dynamical models there is also a small shift in the position of 
this maximum towards lower energy. 



4 Summary and conclusions 

The main objective of this article was to explore the influence of the expan- 
sion dynamics, the equation of state and the nature of deconfinement phase 
transition on strangeness production in heavy ion collisions. 

We have discussed and formulated different models for a phase transi- 
tion in hight density QCD matter. The thermodynamical properties of these 
models and the role of the order of the phase transition as well as the in- 
teractions between the particles has been analyzed. We have addressed the 
question of the Gibbs construction of the phase transition in the presence of 
two conserved charges and emphasized the problem of causality and thermo- 
dynamical consistency. 

The strangeness separation in the transition region from the quark-gluon 
plasma to the hadronic phase was also studied. The asymmetry in the relative 
concentration of strange and anti-strange quarks in the hadronic and quark- 
gluon component in the phase coexistence region was found in all models that 
exhibit a phase transition. However, the largest effect was observed in the 
mixed-phase model with a crossover-type deconfinement phase transition. 

The differences in equilibrium thermodynamics of the models were stud- 
ied on the dynamical level. We have shown that the hydrodynamical expan- 
sion trajectories of the fireball in the T-fi B plane are very sensitive to the 
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equation of state. We considered the effect of the different expansion paths 
on strangeness production. Our detailed analysis show that there is almost 
no sensitivity of strangeness observables on the equation of state or on the 
expansion trajectories. This was demonstrated for several strange particle 
excitation functions. 

To relate the model predictions with experimental data we have extended 
our study to a canonical formulation of strangeness conservation. We have 
discussed the phenomenological limitations of our dynamical models and the 
possible extension needed to provide a quantitative description of the ob- 
served particle yields in heavy ion collisions. 

Exact strangeness conservation substantially reduces the strange particle 
yields in heavy ion collisions for < 10 GeV. For higher energies a mod- 
erate suppression is also found if the beam-energy dependence of the volume 
parameter V c is taken into account. We have shown that the assumption 
that V c is the volume of the initially produced Lorentz contracted fireball 
may lead to a negative slope in the energy dependence of the K + /ir + ra- 
tio. However, within considered models, the almost singular behavior of the 
excitation function near Ei^ < 20 GeV for K + /tt + ratio found recently by 
the NA49 collaboration jJE] was not reproduced. Simplified hydrodynamics 
with the assumption of a shock-like particles freeze-out in heavy-ion col- 
lisions results in a very smooth behavior of the strange particle excitation 
functions. 
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